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We report joint Lagrangian velocity and temperature measurements in turbulent thermal 
convection. Measurements are performed using an improved version (extended auton¬ 
omy) of the neutrally-buoyant instrumented particle (Shew et al. ]2007t that was used by 
[Gasteuil et al. (20071 to performed experiments in a parallelepipedic Rayleigh-Benard 
cell. The temperature signal is obtained from a RF-transmitter. Simultaneously, we 
determine particle’s position and velocity with one camera, which grants access to the 
Lagrangian heat flux. Due to the extended autonomy of the present particle, we obtain 
well converged temperature and velocity statistics, as well as pseudo-eulerian maps of 
velocity and heat flux. Present experimental results have also been compared with the 
results obtained by a corresponding campaign of Direct Numerical Simulations and La¬ 
grangian Tracking of massless tracers. The comparison between experimental and numer¬ 
ical results show the accuracy and reliability of our experimental measurements. Finally, 
the analysis of lagrangian velocity and temperature frequency spectra is shown and dis¬ 
cussed. In particular, we observe that temperature spectra exhibit an anomalous 
frequency scaling, likely representing the ubiquitous passive and active scalar behavior 
of temperature. 


Key words: You can use up to three keywords, but you have to choose them from 
a list of possible alternatives when you submit the paper Rayleigh-Benard convection, 
large-scale circulation, smart particle, turbulent convection, Lagrangian measurements. 


1. Introduction 

Thermal convection occurs in many industrial and geophysical applications, ranging 
from heat exchangers or nuclear/chemical reactors to atmspheric circulation. In most 
of these cases, the flow is highly turbulent and transports heat very efficiently. Never¬ 
theless, understanding and modeling the local and global properties of the temperature, 
velocity and heat-flux helds in these situations is still a challenge. To analyze turbu¬ 
lent thermal convection in a laboratory, we choose the Rayleigh-Benard configuration: 
a horizontal layer of fluid confined between a cooling plate above and a heating plate 
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below. In this flow configuration, the temperature gradient is confined in the thermal 
boundary layers close to the heating and cooling plates. Due to the strong mixing, 
a nearly homogeneous temperature distribution is observed in the bulk. The 
driving force of the flow is measured by the Rayleigh number: 

= ( 1 , 1 ) 

where H is the height of the cell, g is the gravitational acceleration, a is the constant 
pressure thermal expansion coefficient, AT = Th — Tc is the difference of temperature 
between the heating and the cooling plate, v is the kinematic viscosity of the fluid, and k 
is its thermal diffusivity. The Prandtl number expresses the ratio between viscous 
and thermal dissipation: 


Pr = 

K 


( 1 . 2 ) 


A further input parameter is the aspect ratio F, i.e. the ratio between the horizontal 
and the vertical size of the cell. The response of the system is represented by the Nusselt 
number, which compares convective and conductive heat flux: 


Nu = 


QH 

A AT’ 


(1.3) 


where Q is the global heat flux, and A is the thermal conductivity of the fluid. 

Assuming locally homogeneous and isotropic turbulence, for sufficiently high Reynolds 
and Peclet numbers, passive scalars and velocity spectra follow the well-known Kolmogorov- 
Obukhov laws ( Kolmogorov|194T Monin & Yaglom pOOTt . However, in thermally driven 
flow, temperature is not a passive scalar and the similarity theory requires a further 
generalization. It has been argued that for small scales a thermally stratihed fluid can 
be considered locally stationary and homogeneous, but not isotropic and axially sym¬ 
metric relative to the vertical direction ( jBolgiano 1959 Obukhov 1959 1 . Within this 
framework, the scalings become dependent on the Bolgiano-Obukhov (B059) length- 
scale measuring the importance of the thermal stratification. 


Ln = e', 




(1.4) 


where ct and e„ are the temperature and kinetic energy dissipation rate, respectively. Lg 
characterizes the minimum length scale of inhomogeneities beyond which stratification 
should be taken into account. The ordinary Kolmogorov spectrum scalings (fc“®/^) are 
expected to be recovered at scales smaller than the Bolgiano length. On the other hand, 
\i Lg becomes much larger than the external turbulence length scale Tq; the effect of 
the mean flow becomes important and the similarity theory does not apply. Yet, all the 
theoretical predictions and scalings are assessed for regions of space far enough from 
solid boundaries. In general, stratification affects a certain range of scales, which cannot 
be considered locally isotropic. Corrections to the velocity and temperature correlations 
and spectra may be universal and may be in principle determined empirically. For strat¬ 


ified flows, Bolgiano (1959) established a theoretical framework to determine 
the asymptotic form of these functions for scales much greater than Lg. These 
Bolgiano-Obukhov (B059) scalings predict spectra that are steeper for the velocity and 
milder for the temperature compared to those given by Kolmogorov-Obukhov scalings. 
Nevertheless, the ultimate picture on velocity and temperature scalings in thermal con¬ 


vection is far from being obtained (Lohse & Xia 20101. Experimental measurements at 
Ra « 10^° and Pr = 4.4 exhibit a K41 behaviour for velocity (Zhou et al.||2008) and for 

























Temperature and velocity Lagrangian measurements in turbulent convection 3 


velocity and temperature ( Sun et aL|2006 ) structure functions. However, direct numerical 
simulations (DNS) of thermal convection on a cilindrical domain ( [Kunnen et al. 2008) 
show a B059 scaling for temperature structure functions in the radial direction and for 
vertical velocity structure functions in the longitudinal direction (Ra = 10®, Pr = 1). 

Even though remarkable progresses have been made towards a deeper understanding 


of scaling laws between control and response parameters (Grossmann & Lohse 2000 
Chavanne et a/.|2001 Stevens et a/.|2013 |, new investigations are required to clarify some 


of the open issues (Lohse & Xia 2010). Most of literature studies on Rayleigh-Benard 
turbulence have been focused on eulerian measurements of velocity and temperature 
distribution (Xia et al. 2003 Tilgner et al. 1993), with the aim of characterizing the 


behavior of the local heat flux, Nu, as a function of the Rayleigh number, Ra (Shang et al. 

[2004 ). Although the mean velocity in the central region of the convection cell 
is homogeneous and close to zero, the velocity root-mean-square is neither 


homogeneous nor isotropic (Xia et al. 


2003 


Qiu et al. 2000). These flow 


inhomogeneities prevent from using the frozen-flow hypothesis (Taylor 1938): 
this makes the connection between time-domain measurements and space- 
domain predictions difficult. 

More recently, improvements in computing power and storage capacities has allowed 
the appearance of lagrangian studies of turbulence, which naturally provide useful 
informations on transport mechanisms. The Lagrangian description of turbulence 
has significantly contributed to our current comprehension of transfer pro¬ 
cesses dToschi & Bodenschatz||2009|) . In particular, a number of experimente 


et al. 

2001 

et al. 

o 

o 


[2001 |Porta et al. 2001 Voth et al. 2002) and numerical (Yeung 2002[ Biferale 


2004) studies focused on velocity and acceleration statistics in homogeneous and 


isotropic turbulence. The first numerical lagrangian studies of thermal convection were 
those of Schumacher (2008, 2009), which were specifically focused on pair dispersion 


and on acceleration statistics. Although the tracer motion was largely anisotropic due 
to the vertical buoyancy, pair dispersion was close to the homogeneous and isotropic 
turbulence. Acceleration and temperature statistics showed a non-gaussian behaviour 
characterized by a large intermittency (higher in the horizontal directions). Relevant to 
the present work was also the finding of the non-symmetric behavior of the heat trans¬ 
port probability density function. From an experimental point of view, measurements 
of turbulent thermal convection in a Lagrangian framework are relatively scarce. Only 
recently, Ni et al. (2012) used three-dimensional particle tracking velocimetry to analyze 
velocity and acceleration statistics in turbulent Rayleigh-Benard convection. In partic¬ 
ular, they observed a gaussian and a stretched exponential distribution for 
the probability density function of velocity and acceleration in the centre of 
the cell. 

From the above review, phenomenological and statistical analyses of turbulent Rayleigh- 
Benard convection in a Lagrangian frame of reference appear not yet complete and require 
further investigation. This is exactly the purpose of the present study. In this work, we 
present both experimental and numerical measurements of temperature and velocity in 
a Lagrangian frame. For the experiments, we improved the neutrally-buoyant instru¬ 


mented particle presented in Shew et al. (2007) and already used and tested by Gasteuil 


et al. (2007). This smart particle explores a rectangular Rayleigh-Benard cell filled with 


water. We compare Eulerian maps obtained from our Lagrangian data with PIV mea¬ 
surements, to show that the particle samples correctly the entire flow (and to deduce 
pseudo-Eulerian maps of temperature and thermal flux). At the same time, we per¬ 
form Direct Numerical Simulations (DNS) and Lagrangian particle tracking 
of massless tracers in turbulent Rayleigh-Benard convection. These numerical 
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Figure 1. (a) Sketch of the convection cell and of the mobile sensor. The six dark dots in the 
plates show the location of the PTIOO temperature sensors, (b) Photograph of the instrumented 
particle (open). We see the batteries place. 


simulations correspond to the ideal case of a particle of zero diameter and, therefore, they 
can be viewed as a ’’conceptual experiment” built to highlight possible finite-size effect 
of the smart-particle in the experiments. Velocity, temperature and heat flux statistics 
obtained from experiments and numerical simulations are compared and discussed. 


2. Smart particle, experimental setup and numerical simulation 
method 

2.1. Smart Particle 

The mobile sensor consists in a 2.1 cm in diameter capsule containing temperature instru¬ 
mentation, a radio-frequency emitter, and two batteries. Four cylindrical thermistors (0.8 
mm in length, 0.4 mm in diameter, 230 kfl, response time 0.6 s in water) are mounted 
on the capsule wall protruding 0.5 mm into the surrounding flow (see sketch in zoom 
of figure [^a) and photograph figure [^b)). A resistance controlled oscillator is used to 
create a square wave whose frequency depends on the average of the four measured tem¬ 
peratures. This square wave is used directly to modulate the frequency of a radio wave 
generated by the radio-frequency emitter. The temperature signal is recovered on the fly 
by a stationary receiver. The capsule has been redesigned (compared to the one described 


by Shew et al. (20071) in order to have the four thermistors at the equator and a simpler 


handling. A new shell has been conceived to contain two batteries and to extend the 
emission time which can now reach up to 1000 turnover times. At the same time the tra¬ 
jectory of the particle is recorded with a digital camera placed in front of the large face of 
the cell (due to cell dimensions, we assume that the mean flow is quasi bi-dimensional). 
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AT{°C) 

Ra 

Nu 

13.15 

3.5 X 10^° 

230 

18.60 

5.0 X 10^° 

244 

22.90 

6.2 X 10^° 

264 


Table 1. Parameters used for acquisitions in the experiments. 



Figure 2. Sketch of the numerical simulation domain. Temperature contours for Ra = 10® 
are also shown (red indicates regions of high temperature whereas blue indicates regions of low 
temperature). 

2.2. Experimental setup 

Our convection cell is a 10.5 cm-thick 41.5 cm x 41.5 cm rectangular cell with 2.5 cm- 
thick PMMA walls (see sketch figurej^a)). Both plates consist in 4 cm-thick copper plates 
coated with a thin layer of nickel. The bottom plate is Joule-heated while the top plate is 
cooled with a temperature regulated water circulation. Plate temperatures are controlled 
by PT 100 temperature sensors. We work with deionized water. The bulk temperature is 
fixed between 37.05°C and 38.35°C in different experiments. The corresponding Prandtl 
number are in the range 4.62 — 4.49. Main parameters are grouped in table 

2.3. Numerical simulations 

Direct Numerical Simulations are performed to complement our experimental results. We 
consider an incompressible and Newtonian turbulent flow of water confined between two 
rigid boundaries. Streamwise, spanwise and wall-normal coordinates are indicated by x, y 
and z, respectively. The bottom wall is kept at uniform high temperature {Th), whereas 
the top wall is kept at uniform low temperature (T),). The size of the computational 
domain is x Ly x Lz = 4 ,ttH x AttH x 2H (in x, y and z, respectively), where H is the 
half-channel height. A sketch of the computational domain/flow conditions is presented 
in Figj^ 

The imposed temperature difference AT = (Th — Tc) between the bottom and the top 
wall induces an unstable buoyancy effect within the flow field (the acceleration due to 
gravity g acts downward along z). Mass, momentum and energy equations in dimension¬ 
less form and under the Boussinesq approximation are: 


V • u = 0, 


(2.1) 
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| + (u.V)0 = 


--V^e, 


(2.3) 


V PrRa 

where Ui is the component of the velocity vector, 0 is the dimensionless tempera¬ 
ture 9 = [T — Tref)/‘^T, p is pressure, whereas is the buoyancy force (acting in 
the vertical direction only) that drives the flow. Eqs. 2.1|2.3 have been obtained using 
h as reference length, Uref = ^JgaAT/2H as reference velocity, T^ef = (T^ + Tc)/2 as 
reference temperature and p = pgaHAT/2 as reference pressure. Density p, kinematic 
viscosity v, thermal diffusivity k and thermal expansion coefficient a are evaluated at a 
mea n fluid temperature of ~ 30°C. The Prandtl and the Rayleigh numbers in Eqs. [H} 
2.3 are defined as Pr = v/k and Ra = {gaAT{2H)^)/{lyk), respectively. In the present 
study, we keep the Prandtl number Pr = 4 and we vary the Rayleigh number between 
Ra = 10^ and Ra = 10®. Periodic boundary conditions are imposed on velocity and 
temperature along the streamwise x and spanwise y directions; at the walls, no slip con¬ 
ditions are enforced for the momentum equations while constant temperature conditions 
are adopted for the energy equation. The resulting set of equations are discretized using a 
pseudo-spectral method based on transforming the field variables into wavenumber space, 
through Fourier representations for the periodic (homogeneous) directions x and y, and 


Chebychev representation for the wall-normal (non-homogeneous) direction z (see Zonta 


et al. (20121; Zonta & Soldati (2014| for details). We used up to 512 x 512 x 513 grid 
points to discretize the computational domain. We injected Np = 1.28 • 10^ lagrangian 
tracers and we computed their dynamics as 


= u(xp(t),t) dp = 9 {xp{t),t). 


(2.4) 


with Xp the tracers position and 9p their temperature. Velocity and temperature at 
particle position are obtained by 6*^ order Lagrange polynomials. Time advancement for 
the Lagrangian tracers is achieved using a 4‘^ order Runge-Kutta scheme. 

As discussed above, numerical simulations are carried out in a laterally-opened (along x 
and y) domain configuration, which is similar to that considered in the theory (where the 
vertical direction is not homogeneous, because of the buoyancy, while the other two are 
homogeneous and isotropic). Experiments must be performed in a closed cell and some 
anisotropic effects may appear also in the other directions. Moreover, passive tracers in 
numerical computations are free to explore the entire domain, whereas the smart-particle 
cannot access the small region near-to-the-wall, given its dimensions. These differences 
might have some effects, as will be discussed in the following. 


3. Lagrangian measnrements 

To follow accurately the flow, the capsule and fluid density are carefully matched 
within 0.05%. This is the main difficulty of the experiment. On figure we plot the 
skewness of the distributions of the horizontal and vertical positions while varying mean 
temperature of the cell. Due to the flow symmetry, we assume that the particle ideally 
matches the fluid density when both horizontal and vertical skewness are close to 0, 
which happens for a temperature of 37.5°C. If the mean temperature is shifted of few 
tenth of degrees, symmetry is broken while the particle becomes less neutrally-buoyant. 
The effect is more dramatic on the horizontal position. To explain it, we assume a particle 
denser than the fluid. When traveling along the top plate, it is easily advected by a cold 
plumes: the particle goes downwards earlier during the travel. When traveling along the 
bottom plate, it is more difficult for a plume to advect the particle: it goes upwards only 
when reaching the corner of the cell. As a consequence, the average horizontal particle 
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Figure 3. Skewness of vertical and horizontal positions of the smart particle versus the 

temperature of the cell center. 
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Figure 4. (a) Flow temperature along the particle trajectory at Ra = 5.0 X 10^°, six-hour 
acquisition, (b) Corresponding Lagrangian vertical Nusselt Nu^. Trajectories are undersampled 
for visibility purposes. 


trajectory is shifted close to the wall where hot plumes rise. This reasoning holds also 
for a less dense particle and explain the large skewness of the horizontal position. 

We performed measurements from 6 to 20.3 hours. Figure |^a) shows an example of 
a temperature measurement along the particle trajectory. In this case, Ra = 
5.0X10^° while the acquisition time is six hours. Globally, the particle describes a loop in 
the counterclockwise direction (the rotation direction can change for other acquisitions). 
Its mean speed is 1 cm/ s. Close to walls and plates, its speed is typically 2 to 3 cm/ s. 
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The thermal boundary layer thickness can be computed as: 


< 5 « = 


H 

2Nu' 


(3.1) 


which corresponds to a thickness of less than 1 mm for the considered Rayleigh num¬ 
bers, which is not directly accessible to our particle. This is why we do not observe 
large temperature gradients close to the top and bottom plates. However, we 
clearly detect hot and cold jets near the right and left walls, respectively. Temperature 
fluctuations in the cell (outside the boundary layers) are typically 1°C and are likely 
dne to advection by plnmes. Nevertheless, joint measurements of temperature 
and trajectory give indications that the particle movement is more influenced 
by the mean wind than by plumes, except close to the vertical walls - inside hot 
and cold jets. 

Knowing the velocity and temperature of the particle, we can compute the Lagrangian 


thermal flux. We use a normalized vertical Nusselt number (Ching et al. 2004 Grossmann 
fc Lohse|[20Ml ): 


= 1 + ^ im - {T{t))t) vAt), (3.2) 

where T{t) is the instantaneous temperature measured by the particle, {T(t))t is its 
mean along the trajectory and Vz{t) is its vertical velocity. Figure Qb) shows the Nu^ 
along the trajectory. Most of the vertical thermal transfer is symmetrically concentrated 
in the hot and cold jets, corresponding to plumes with T{t) — {T{t))t and Vz{t) having 
the same sign. We also observe that the vertical flux is highly asymmetric towards 
positive value up to 30 times the average vertical heat flux {Nu^)t = 138. 


4. Pseudo-Eulerian maps 

4.1. Methodology and PIV measurements 

The velocity distribution in the central region of the convection cell is close to that of a 
solid body rotation, a situation that hinders the particle from exploring easily this 
region. To obtain a correct resolution of the central region we performed very long exper¬ 
iments. With more than twenty-hour measurement, we have enough data in the whole 
cell (including the central zone) to compute pseudo-Eulerian maps of several quantities. 
To do this, we divide our cell in 1.04 cm x 1.04 cm squares and we compute the 
average of the considered quantity in each cell. The accuracy of the method has been 
evaluated through comparison of the pseudo-Eulerian velocity field (figure]^a)) with 
the results obtained by Particle Image Velocimetry (PIV) measurements in the same 
cell (figure l^b)) at similar Rayleigh numbers. We performed PIV measurements with a 
1.2 W, Nd:YV04 laser. Flow was seeded with Spherical HOPS glass beads of 1.10±0.05 
in density and of 12 pm average diameter. Twelve-hour acquisitions with one picture 
pair (frequency acquisition 20 Hz) every ten seconds were used to compute mean veloc¬ 
ity fields. We used CIVx (Fincham & Delerce 2000) free software for analysis. Several 
passes are applied to picture pairs. For the first one, we cut out pictures in 30x30 pixels^ 
elementary boxes with 50% overlap. Search zones was one and a half larger. 

The velocity module obtained from PIV is in good agreement with observations from 


Xia et al. (2003) for a similar cell. The velocity measured with the smart particle is 


slightly lower than that obtained with the PIV technique, probably due to several 
effects: a small difference in the Rayleigh number, particle size and inertia, and fluid 
slip on the particle surface. 
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Figure 5. (a) Pseudo-Eulerian module velocity field from Lagrangian data, Ra = 5.0 x 10^°. 
(b) Module velocity field obtained by PIV, Ra = 5.6 x 10^°. 



4.2. Temperature and thermal flux maps 

Our Lagrangian method can be efficiently used to obtain pseudo-Eulerian map of tem¬ 
perature in the whole cell. To our knowledge, this is the first Lagrangian experiment 
giving the whole temperature map for this range of Rayleigh and Prandtl numbers. The 
pseudo-Eulerian temperature field is plotted in figure [^a) . We observe that the flow is 
homogeneous in the bulk, whereas hot and cold jets dominate the regions close to the 
walls. Deviations of temperature from the bulk value are in particular seen in the top-left 
and bottom-right corners and are likely due to the effect of buoyant plumes driving the 
smart particle along the vertical walls. As suggested by Scagliarini et al. (2014), 
the mean wind acts stabilizing the boundary layers and reducing the plume 
emission activity. Moreover, the intense mixing makes the plume temperature close 
to the bulk temperature along the plates but in the top-left and bottom-right corners 
(where the mean wind is blocked and mixing cannot be observed). This situation is 
well represented in figure [^b), where we observe large fluctuations in the 
corners but fewer fluctuations in the vertical jets and close to the plates. The 
observed slight asymmetry can be explained by a non-perfect particle-to-fiuid 
density matching. With velocity and temperature joint measurements, this is the only 
Lagrangian experimental technique that gives a pseudo-Eulerian thermal flux map. We 
observed that Nu^ is large inside both cold and hot jets, where vertical ve¬ 
locity Vz is large, and indicates that the spatial distribution of Nuf is chiefly 
influenced by plumes. 


4.3. Velocity fluctuations 

Now we discuss velocity root mean square (RMS) maps obtained by our smart parti¬ 
cle in connection with the results obtained by PIV measurements. From the velocity 
pseudo-Eulerian maps, we can compute smooth velocity mean fields by interpolation. 
This velocity value at each point {x,z) is indicated as vf{x,z) where i = x, z is the 
horizontal or vertical velocity, respectively. Thus, in each square s described above, for 
all {xs,ys) Lagrangian coordinates of the particle trajectory included in the square we 
have: 
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Figure 7. Vertical Nusselt Nu^ pseudo-Eulerian map at Ra = 5.0 x 10^°. 



RMS 

^i,s 



(4.1) 


where v("(xs, Zg) represents the Lagrangian velocity events inside the considered square 
and {.)s is the average of these events. This is similar to the velocity RMS computed from 
PIV measurements for all {x, y) in the cell: 


vp^^{x,z) = 


Vi{x,Z,t) - {Vi{x,z,t))tf)t. 


(4.2) 


In these two definitions, both vf{x, z) and {vi{x, z, t))t represent the mean flow. Figure 
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Figure 8. Top-left: vertical velocity RMS pseudo-Eulerian map; top-right: horizontal velocity 
RMS pseudo-eulerian map. Bottom-left: vertical velocity RMS obtained by PIV; bottom-right: 
horizontal velocity RMS obtained by PIV. Ra — 5.0 x 10^° for pseudo-Eulerian maps and 
Ra = 5.6 X 10 “ for PIV fields. 


1^ compares vertical (left) and horizontal (right) velocity fluctuations from the smart 
particle (top) and from PIV (bottom). First, we note that RMS values recorded by the 
particle are slightly smaller compared to those evaluated by PIV. This is mainly due 
to the filtering effect played by the particle on small scales fluctuations. There is also 
a secondary effect due to a small difference in the value of Ra {Ra = 5 • 10^° for 
the Lagrangian particle and Ra — 5.6 • 10^° for PIV measurements). The two regions 
characterized by the largest fluctuations are those where vertical plumes impinge on the 
horizontal walls and induce large turbulence patches. In our Lagrangian measurements, 
we also observe significant fluctuations along the horizontal plates: these fluctuations, 
which are not visible in the PIV measurements, are due to particle rebounds on the 
horizontal walls. When the particle hits the wall, a positive vertical velocity becomes a 
negative vertical velocity and vice versa. Differently, when the fluid approaches the wall, 
a large part of its vertical velocity is converted into horizontal one. Consequently, for 
each rebound, we can suppose that the contribution of the particle to the RMS is twice 
that of the fluid. This leads to a 50% increase of the RMS value, which corresponds 
to the increase observed between pseudo-Eulerian and Eulerian vertical velocity RMS 
maps. 
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5. Flow statistics 

5.1. Probability density functions 

In figure l^a) we show the probability density function (PDF) of the temperature fluctu¬ 
ations recorded by our instrumented particle at Ra = 5 • 10^° for different measurement 
time {t = 2,4,8,16,20.5 hr). We note that the shape of the PDF becomes increasingly 
smooth for increasing measurement times. Present results based on 20.5-hour measure¬ 


ments substantially improve previous results of Gasteuil et al. { 2007) obtained with a 


2-hour acquisition. Nevertheless, the global shape is conserved, and confirms the 
overall quality of the results by Gasteuil et al. ( }2007| ). The symmetric shape of 
the PDF is a further signature of the quality of the results (good buoyancy neutrality). 
The PDF has sharp exponential tails, in agreement with previous Eulerian measure¬ 


ments (Belmonte et al. 19941 performed far from the boundary layers and using air at 
Ra = 4.8 X 10 b However, tails are wider due to the passage of the particle in hot and 
cold jets. This suggests that Eulerian measurements in Rayleigh-Benard convection are 
delicate: since the flow is largely inhomogeneous, the position at which the measurement 
is taken is fundamental. The generalization of the behaviour of the entire cell based on a 
local measurement requires a lot of care. In Fig. |^b), we compare experimental results 
of the PDF of temperature fluctuations obtained in the 20.5-hour measurement with our 
results from numerical simulations at different Rayleigh numbers. Results are normalized 
by the corresponding temperature standard deviation, std(T). Even though numerical 
experiments are carried out at smaller Rayleigh number, the agreement between experi¬ 
ments and simulations is satisfactory when Ra > 10®. Deviations between experimental 
and numerical results are observed only for (T— {T{t))t)/std{T) > 3, hence highlighting 
the accuracy of the present experiments. Indeed, our numerical simulations can be seen as 
an ideal experiment, since we sample the flow-field with massless (pointwise) fluid tracers 
(no size/inertia effect of the Lagrangian probes). Within this framework, the difference 
between experiments and simulations can give indications on the effect of the size of our 
smart particle: due to the finite-size, our smart particle acts as a filter for small/short 
space/time scale events, and numerical experiments show larger tails. 

As previously mentioned, with our instrumented particle we are able to record simul¬ 
taneously velocity and tempe ratu re, hence we can compute the local value of the vertical 
Nusselt number Nu^. Figure [l0l(a) shows the histogram of the vertical Nusselt for three 
different Rayleigh numbers. The most probable value is Nujt = 0, corresponding to the 
time during which the particle has a horizontal trajectory or is advected by the mean 
wind far from the walls (hence far from plumes), where no significant vertical heat flux is 
observed. Interestingly, there is a larger positive tail, which is the signature of near wall 


intense events (Gasteuil et al. 2007). When Ra increases, the shape of the histogram 
does not change, but we observe more intense events at higher Ra. In figure |10[ b) we 
compare our experimental and numerical results of PDF{Nu^). Results are normalized 
by their corresponding standard deviation, std{Nu^). We observe a very good agreement 
between experiments and simulations over the entire range of measured Nu^/std{Nu^). 
Deviations are seen only for extreme and rare events, Nujt/std{Nu^) > 7, due to the 
filtering effect of the smart particle size on velocity/temperature fluctuations. 


5.2. Spectral analysis 

We conclude our discussion with a spectral analysis of velocity and temperature La¬ 
grangian time series recorded at different Rayleigh numbers. Figures [Tl|(a) , [TT] (b) and 
12 show the frequency spectra (for six-hour acquisitions) of vertical velocity, horizontal 
velocity and temperature, respectively. We observe three distinctive characteristics. First, 
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Figure 9. (a) Probability Density Function (PDF) of the temperature fluctuations T — (T{t))t 
at Ra — 5.0 x 10^'^ for different measurement times, (b) Comparison between experimental 
and numerical results of the PDF of temperature fluctuations, T — {T{t)}t, normalized by the 
corresponding standard deviation. Solid line represents experimental data whereas dashed lines 
correspond to DNS data at different Ra. 


(a) 



Nu^ 



Nu^/std{Nu^) 


Figure 10. (a) PDF of the vertical thermal flux Nu^ for different Rayleigh numbers, six-hour 
measurements, (b) Comparison of the vertical Nusselt PDF Nu^ normalized by the standard 
deviation. Solid line represents experimental data whereas dashed lines correspond to DNS ones. 


a peak appears at a frequency / ~ 1.25 x 10“^ Hz (80 s): it is consistent with the typical 
large scale turnover time. Yet, when Ra increases, the spectrum shifts upwards, indi¬ 
cating a more intense dynamics. Finally, a cut-off occurs at 0.15 Hz for velocity and 0.5 
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(b) 



Frequency [Hz] 


Figure 11. (a) Vertical and (b) horizontal velocity spectra at different Rayleigh numbers. 


Hz for temperature. This latter observation indicates a finite-size particle effect, which 
is less important for temperature owing to the thermistors lower size. 

From the proposed scaling laws (lines in figures inE3 we note some interesting fea¬ 
tures. A f~^ power law characterizes the vertical velocity, whereas a power law 

characterizes the horizontal velocity. A similar flow anisotropy has been also observed 
by Qureshi et al. (2007) in a Von-Karman flow. From the above observations, we can in¬ 
fer the following physical interpretation. The particle has a large vertical velocity when it 
enters hot and cold jets along the lateral walls, whereas it has a large horizontal velocity 
when it is driven by the mean wind along the top and bottom walls. Since the hot and 
cold jets are characterized by more intense fluctuations, the corresponding spectra have 
a weaker slope. After the cut-off, a classical f~^ power law appears for both spectra. 
Concerning the temperature spectra, we observe a power law. Note that due to the 
anisotropic and inhomogeneous flow condition we do not find a f~^ slope. Yet, temper¬ 
ature can not be considered as a pure passive scalar, in particular inside hot and cold 
jets. Unfortunately, at present we do not have a quantitative explanation for this power 
law. Results from numerical simulations (not shown here) differ from the present exper¬ 
imental results and essentially follow the classical K41 scaling. This difference is mainly 
due to the different flow configuration: our simulations are run in a domain with top and 
bottom walls but no side walls, whereas experiments are run in a square enclosure. 


6. Conclusion 

In this work, we have used an improved version of our smart particle to perform 
new experimental measurements in Rayleigh-Benard convection at different Rayleigh 
numbers (up to Ra = 6.2 • 10^*^). To corroborate the experiments and isolate possible 
finite-size effects, we have also carried out Direct Numerical Simulations of Rayleigh- 
Benard turbulence with Lagrangian tracking of massless tracers. 

In the experiments, due to the extended autonomy of the particle, we are able to sample 
velocity and temperature at the particle position for long periods of time, up to twenty 
hours. This long data recording allows not only to acquire long time-series describing the 
temporal evolution of a turbulent flow, but also to build pseudo-Eulerian maps of the 
flow field and to compute converged PDF and spectra. The particle trajectory is driven 
by the interaction between the mean large scale circulation (along the horizontal walls) 
and the thermal plumes generating vertical hot and cold jets. Velocity and temperature 
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Frequency [Hz] 

Figure 12. Temperature spectra at different Rayleigh numbers. 


fluctuations are essentially concentrated in the bottom-right and top-left corners, and 
denotes strong turbulence events in these regions. 

Upon comparison between experimental and numerical results, we are able to demon¬ 
strate the accuracy of our experimental technique in recovering all the fundamental 
statistical features of the flow. We finally computed velocity and temperature frequency 
spectra. Interestingly, horizontal and vertical velocity spectra exhibit different scaling 
(/-^ and respectively), as a consequence of the strong flow anisotropy in the ver¬ 
tical and horizontal directions. For temperature, we observe a steep scaling, likely 

representing the hybrid passive and active scalar behavior of temperature. 

A further study of these observations will be the subject of a future paper. As the 
flow is dominated by a mean vortical structure, we can remove it from the recorded 
signal to study the behavior of turbulence fluctuations and the interactions between 
fluctuations and mean structure along the line proposed by Machicoane & Volk (2014) 
for Von-Karman flows. 
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